High-Precision Lattice QCD Confronts Experiment 
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We argue that high-precision lattice QCD is now possible, for the first time, because of a new 
improved staggered quark discretization. We compare a wide variety of nonperturbative calculations 
in QCD with experiment, and find agreement to within statistical and systematic errors of 3% or less. 
We also present a new determination of a^^{AIz); we obtain 0.121(3). We discuss the implications 
of this breakthrough for phenomenology and, in particular, for heavy-quark physics. 



PACS numbers: 11.15.Ha,12.38.Aw,12.38.Gc 

For almost thirty years precise numerical studies of 
nonperturbative QCD, formulated on a space-time lat- 
tice, have been stymied by our inability to include the 
effects of realistic quark vacuum polarization. In this 
paper we present detailed evidence of a breakthrough 
that may now permit a wide variety of nonperturbative 
QCD calculations including, for example, high-precision 
B and D meson decay constants, mixing amplitudes, and 
semi-leptonic form factors — all quantities of great im- 
portance in current experimental work on heavy-quark 
physics. The breakthrough comes from a new dis- 
cretization for light quarks: Svmanzik-improved stag- 
gered quarks 011111111111113. 

Quark vacuum polarization is by far the most expen- 
sive ingredient in a QCD simulation. It is particularly dif- 
ficult to simulate with small quark masses, such as u and 
d masses. Consequently, most lattice QCD (LQCD) sim- 
ulations in the past have either omitted quark vacuum 
polarization ("quenched QCD"), or they have included 
effects for only u and d quarks, with masses 10-20 times 
larger than the correct values. This results in uncon- 



trolled systematic errors that can be as large as 30%. The 
Symanzik-improved staggered-quark formalism is among 
the most accurate discretizations, and it is 50-1000 times 
more efficient in simulations than current alternatives of 
comparable accuracy. Consequently realistic simulations 
are possible now, with all three flavors of light quark. An 
exact chiral symmetry of the formalism permits efficient 
simulations with small quark masses. The smallest u and 
d masses we use are still three times too large, but they 
are now small enough that chiral perturbation theory is 
a reliable tool for extrapolating to the correct masses. 

In this paper we demonstrate that LQCD simulations, 
with this new light-quark discretization, can deliver non- 
perturbative results that are accurate to within a few per- 
cent. We do this by comparing LQCD results with exper- 
imental measurements. In making this comparison, we 
restrict ourselves to quantities that are accurately mea- 
sured (< 1% errors), and that can be simulated reliably 
with existing techniques. The latter restriction excludes 
unstable hadrons and multihadron states (e.g., in non- 
Icptonic decays); both of these are strongly affected by 
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the finite volume of our lattice (2.5 fm across). Unstable 
hadrons, like the p and the 4>, are constantly fluctuating 
into on-shell or nearly on-shell decay products that can 
easily propagate to the boundaries of the lattice; simi- 
lar problems afflict multihadron states. Consequently we 
focus here on hadrons that are at least 100 MeV below 
decay threshold or have negligible widths (J/i/', TT- ■ • ); 
and we restrict our attention to hadronic masses, and to 
hadronic matrix elements that have at most one hadron 
in the initial and final states. These are the "gold-plated" 
calculations of LQCD — calculations that must work if 
LQCD is to be trusted at all. 

Unambiguous tests of LQCD are particularly impor- 
tant with staggered quarks. These discretizations have 
the unusual property that a single quark field cre- 
ates four equivalent species or "tastes" of quark. "Taste" 
is used to distinguish this property, a lattice artifact, 
from true quark flavor. A quark vacuum polarization 
loop in such formalisms contributes four times what it 
should. To remove the duplication, the quark determi- 
nant in the path integral is replaced by its fourth root. 
This construction introduces nonlocalities that are poten- 
tially worrisome, but much is known about the formalism 
that is reassuring: for example, 

• perturbation theory, which governs the theory's 
short-distance behavior, is correct to all orders; 

• phenomena, such as tt" — > 27, connected with chi- 
ral anomalies are correctly handled (because the 
relevant (taste-singlet) currents are only approxi- 
mately conserved); 

• the CP violating phase transition that occurs when 
•niu + TTLd < does not occur in this formalism, but 
the real world is neither in this phase nor near it; 

• the nonperturbative quark loop structure is cor- 
rect up to short-distance taste-changing interac- 
tions, which are perturbative; these interactions 
are suppressed by a'^ag and can be systematically 
removed Q; or they can be removed after the 
simulation using modified chiral perturbation the- 
ory (laUll]. 

To press further requires nonperturbative studies. The 
tests we present here are among the most stringent non- 
perturbative tests ever of a staggered quark formalism 
(and indeed of LQCD). 

The gluon configurations that we used, together with 
the raw simulation data for pions and kaons, were pro- 
duced by the MILC collaboration; heavy-quark propaga- 
tors came from the HPQCD collaboration. The lattices 
have lattice spacings of approximately a — 1/8 fm and 
a — 1/11 fm. The simulations employed an 0{a^) im- 
proved staggered-quark discretization of the light-quark 
action '4], a "tadpole-improved" 0(0^ as) accurate dis- 
cretization of the gluon action 0{a'^,v^) improved 
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FIG. 1: LQCD results divided by experimental results for 
nine different quantities, without and with quark vacuum po- 
larization (left and right panels, respectively). The top three 
results are from our a = 1/11 and 1/8 fm simulations; all 
others are from a = 1/8 fm simulations. 



lattice version of NRQCD for b quarks 0| , and the Fer- 
milab action for c quarks |l4| . Several valence u/d quark 
masses, ranging from ms/2 to rris/S, were needed for ac- 
curate extrapolations, as were sea u/d masses ranging 
between TOs/2 and iris/Q. Only u, d and s quark vacuum 
polarization was included; effects from c, 6 and t quarks 
are negligible (< 1%) here. 

To test LQCD, we first tuned its five parameters to 
make the simulation reproduce experiment for five well- 
measured quantities. The five parameters are the bare 
u and d quark masses, which we set equal, the bare s, 
c and b masses, and the bare QCD coupling. There are 
no further free parameters once these are tuned. 

Setting m„ = simplifies our analysis, and has a 
negligible effect (< 1%) on isospin-averaged quantities. 
We tuned the u/d, s, c, and b masses to reproduce mea- 
sured values of to^, 2m\^ — rn^, fno^, and mx, respec- 
tively. In each case the experimental quantity is approxi- 
mately proportional to the corresponding parameter, and 
approximately independent of the other parameters. 

Rather than tune the bare coupling, one normally sets 
the coupling in LQCD to a particular value, and deter- 
mines the lattice spacing a in its place (after the simu- 
lation). We adjusted the lattice spacing to make the T- 
T' mass difference agree with experiment. We chose this 
mass difference since it is almost independent of all quark 
masses, including, in fact, the 6 mass We could 

equally well have chosen, instead, any of the nine test 
quantities discussed below, with similar results. 

Having tuned all free parameters in the simulation, 
we then computed a variety of experimentally accessible 
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FIG. 2: Chiral fits to LQCD determinations of and fn 
(in GeV) for different values of the valence ii/d-quark mass 
at a = 1/11 fm. 



quantities (in addition to the five used for tuning). Our 
results are summarized in Fig. ^ where we plot the ratio 
of LQCD results to experimental results for nine quanti- 
ties: TT and K decay constants, a baryon mass splitting, 
a Bs~T splitting, and mass differences between various 
J/ip and T states. On the left we show ratios from QCD 
simulations without quark vacuum polarization (n / = 0) . 
These results deviate from experiment by as much as 10- 
15%; the deviations can be made as large as 20-30% by 
tuning QCD's input parameters against different physi- 
cal quantities. The right panel shows results from QCD 
simulations that include realistic vacuum polarization. 
These nine results agree with experiment to within sys- 
tematic and statistical errors of 3% or less — with no free 
parameters. 

The dominant uncertainty in the light-quark quanti- 
ties in this plot (the top four) comes from extrapolations 
in the sea and valence light-quark masses. We used par- 
tially quenched chiral perturbation theory to extrapolate 
pion and kaon masses, and the weak decay constants 
and /if. Chiral perturbation theory was unnecessary 
for correcting the s-quark mass; simple linear interpola- 
tion is adequate, and preferable since chiral perturbation 
theory converges slowly for masses as large as rrig. We 
also kept u/d masses smaller than ms/2 in our fits, so 
that low-order chiral perturbation theory was sufficient. 
Our chiral expansions included the full first-order con- 
tribution and also approximate second-order terms, 
which are essential given our quark masses. We corrected 
for errors caused by the finite volume of our lattice (1% 
errors or less) , and by the finite lattice spacing (2-3% er- 
rors). The former corrections were determined from chi- 
ral perturbation theory; the latter by comparing results 
from the coarse and fine lattices. Residual discretization 
errors, due to nonanalytic taste- violations were 
estimated as 1.9% for /^r and 1% for Jk- Perturbative 



matching was unnecessary for the decay constants since 
they were extracted from partially conserved currents. 
Our final results for fx and agree with experiment to 
within systematic and statistical uncertainties of 2.8%. 
For the n/ = case we analyzed only a = 1/8 fm, but 
extrapolated to the continuum in an approximate way 
based upon our Uf — 3 analysis. 

Fig. El which shows our fits for and Jk, demon- 
strates that the u/d masses currently accessible with im- 
proved staggered quarks are small enough for reliable 
and accurate chiral extrapolations, at least for pions and 
kaons. The valence and sea s-quark masses were 14% 
too high in these simulations; and the sea u/d masses 
were ms/2. 3 and ms/4.5 for the top and bottom results 
in each pair. The dashed lines show the fit function with 
corrected s and sea-quark masses; these lines extrapolate 
to the final fit results. The extrapolations are not large — 
only 4-9%. Indeed the masses are sufficiently small that 
simple linear extrapolations give the same results as our 
fits, within few percent errors. These decay constants 
represent the current state of an ongoing project; a more 
thorough analysis will be published soon. 

The other quantities in the ratio plot, Fig.^ are much 
less sensitive to the valence u/d mass and soft-pion ef- 
fects. Consequently, they are more stringent tests of 
LQCD. The combinations ^M^-Mn and 2Mb, -Mr de- 
pend upon the valence s mass, but the s masses we used 
in our simulations are off by only 10-20% and easily cor- 
rected. The b's rest mass cancels in 2Mb, — Mr, making 
this a particularly clean and sensitive test. The same 
is true of all the T splittings, and our simulations con- 
firm that these are also independent (< 1-2%) of the sea 
quark masses for our smallest masses. The T(P) masses 
are averages over the known spin states; the T{1D) is the 
1^D2 state recently discovered by CLEO 17] . 

It is important to appreciate that our heavy-quark re- 
sults come directly from the QCD path integral, with 
only bare masses and a coupling as inputs — five num- 
bers. Furthermore, unlike in quark models or HQET, 
T physics in LQCD is inextricably linked to B physics, 
through the 6-quark action. Our results strongly suggest 
that effective field theories, like NRQCD, are reliable and 
accurate tools for analyzing heavy-quark dynamics. 

Another important ingredient in high-precision LQCD 
is perturbation theory, which connects lattice results to 
the continuum. We tested perturbation theory by ex- 
tracting values of the coupling from our simulations and 
comparing them with non-LQCD results. We determined 
the renormalized coupling, ay (6. 3 GeV), by comparing 
2"'^-order perturbation theory for the expectation value 
of a 1 X 1 Wilson loop with (exact) values from the simu- 
lations 0, 01 . Results for several sea-quark masses are 
shown in Table ^ the masses become more realistic as 
one moves down the table. 

The QCD coupling is sensitive to the tuning of the lat- 
tice spacing, since this in effect tunes the bare coupling. 
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TABLE I: The QCD coupling av(6.3GeV) from 1 x 1 Wil- 
son loops in simulations with different u/d and s sea-quark 
masses (in units of the physical s mass), and using two dif- 
ferent tunings for the lattice spacing. The first error shown is 
statistical, and the second is truncation error (0(0:^)). 
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We show results for two different tunings: one using the 
T{1P - IS") splitting, and the other using T{2S - IS). 
The two tunings give couplings that are ten standard 
deviations apart and 25% smaller than the physical cou- 
pling when the sea-quark masses are infinite. 

With smaller, more realistic sea-quark masses, the two 
tunings agree to within 1%, and the coupling becomes 
mass independent. Our results, converted to MS and 
evolved perturbatively to scale Mz, imply 

c^^^{Mz) = 0.121 {3), (1) 

which a gree s well with the current world average of 
0.117 (2) []J|. Ours is the first determination from lattice 
QCD simulations with realistic quark vacuum polariza- 
tion, the first with O(a^) improved actions, and the first 
that is verified by a wide range of heavy-quark and light- 
quark calculations; and it is by far the most thorough 
study of the light-quark mass dependence (or indepen- 
dence) of lattice QCD determinations. A more detailed 
discussion will be presented elsewhere. 

The results presented here suggest that we now have 
a reasonably generic and accurate tool for solving a real- 
life, strongly coupled, quantum field theory — for the 
first time in the history of particle physics. Much is re- 
quired to complete the argument. Chiral extrapolations 
for non-strange baryons, for example, are expected to 
be much larger than for pions and kaons, as are finite- 
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FIG. 3: Gold-plated LQCD processes that bear on CKM ma- 
trix elements. eK is another gold-plated quantity. 



volume errors; computations with these hadrons are not 
yet under control. Also a wider variety of tests is im- 
portant. Heavy-quark mixing amplitudes, and semilep- 
tonic decay form factors, for example, are essential to 
the high-precision experiments at B factories; our lattice 
techniques for these require independent tests. The new 
CLEO-c program will be particularly useful for this. 

The larger challenge facing LQCD is to exploit these 
new techniques in the discovery of new physics. Again, 
B and D physics offer extraordinary opportunities for 
new physics from LQCD. There are, for example, gold- 
plated lattice quantities for every CKM matrix element 
except Vtb (Fig. 13). An immediate challenge is to pre- 
dict the D I Dg leptonic and semi-leptonic decays rates to 
within a few percent before CLEO-c measures them. 
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